Giant Rabi splitting in metallic cluster - cavity system 
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We investigate theoretically the photo absorbtion of the cluster of alkali atoms embedded into a 
single mode quantum microcavity. We show that when the energy of the giant plasmonic resonance 
lies close to the energy of the cavity mode, the strong coupling between plasmon and cavity photon 
can occur, which is characterized by mode anticrossing and observation of the doublet structure in 
ptotoabsorbtion. The characteristic values of the Rabi splitting are expected to be several orders of 
magnitude larger then those observed in single quantum dot- cavity systems. 



Introduction Strong coupling between light and mat- 
ter excitations attract the growing interest of the physical 
community. The problem is important not only because 
of the fundamental aspects brought forward by the in- 
teraction of material systems with photons [l[, but also 
due to the possibility of use of the strong coupling phe- 
nomena for the creation of optoelectronic devices of new 
generation, such as polariton lasersl2|, optical logic gates 
[3[, all-optical integrated circuits [4] and others. 

Coupling of zero dimensional (OD) quantum system to 
a single photon mode, which forms the subject of cav- 
ity Quantum Electrodynamics (cQED), is of particular 
importance from this point of view, because of potential 
application of cQED to quantum information processing 
[SKZl. In the domain of the condensed matter, the system 
which recently attracted particular attention consists of 
a Quantum Dot (QD) coupled to a single microcavity 
mode (8l-[l2|. The material excitations in the QD are 
excitons, that is, bound electron-hole pairs. Owing to 
their spatial confinement and energy level discretization, 
they can be brought in strong coupling with the single 
mode of a microcavity, such as that offered by a pillar 
(etched planar cavity) [9|], the defect of a photonic crys- 
tal 10 1 or the whispering gallery mode of a microdisk 
[m . 121, among others. In Refs. (9l-[ll|. such structures 
have demonstrated Rabi doublet in their optical spectra, 
which characterises the mode anticrossing that marks the 
overcome of dissipation by the coherent exciton-photon 
interaction. The characteristic values of Rabi splitting 
for QDs are typically of the order of 10-100 /ieV, which 
is 3-4 orders of magnitude smaller that for the planar 
microcavities with quantum wells. The small value of 
the Rabi splitting makes the achievement of the regime 
of strong coupling in OD systems technically complicated 
task and limits the possibilities of its practical implemen- 
tations. 

The natural question arises: can strong coupling 
regime be observed for material excitations other then 
excitons? One of the natural candidates is collective plas- 
monic excitation in metals. For planar metallic struc- 
tures the dispersion of the surface plasmon lies outside 
the light cone, which makes its direct optical excitation 
impossible [13| and rules out the possibility of the ob- 
servation of any strong coupling effects. This is not 
true, however, for more complicated structures contain- 



ing metallic nanowire arrays [1^ and metallic nanorods 
[15|, for which the effects of strong coupling were shown 
to be extremely pronounced and experimentally observed 
Rabi splitting can be as big as 250 fieV Q. This num- 
ber exceeds by the order of magnitude the values of Rabi 
splitting for exciton-photon coupling in planar inorganic 
microcavities [l6| and is comparable to Rabi splittings 
observed in organic structures [l7| . 

In the present work we consider another type of hy- 
brid metal-dielectric structure, consisting of the single 
metallic cluster embedded inside a single mode photonic 
cavity. The collective motion of electrons in the cluster 
against the positively charged ionic background leads to 
formation of a surface plasmonic mode responsible for the 
appearance of a giant resonance in the photoabsorbtion 
spectra 18-22[. If the energy of the plasmon is close to 



that of the photonic mode of the cavity, the processes of 
multiple resonant emissions and absorptions of photons 
by the cluster can take place and hybrid plasmon-photon 
modes be formed. The goal of the present letter is to 
analyze their influence on photoabsorbtion spectrum of 
the system and to estimate the corresponding values of 
the Rabi splitting. 

The model. To illustrate the onset of the strong cou- 
pling regime in cluster- cavity system we will focus on 
the metallic clusters formed by alkali atoms (Li, Na, K), 
as they are the most exhaustively studied from both the- 
oretical and experimental points of view. The geometry 
of the system we consider is shown on Fig(T] 

The theoretical basis of studying of the collective phe- 
nomena in clusters is provided by the jellium model, in 
which N electrons (one per each alkali atom) are moving 
in the spherically symmetric electrostatic potential Vb{r) 
formed by the uniform background of the positive ions, 
calculated as convolution of the Coulomb electrostatic 
potential with smooth charge distribution. 
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This assumption is generally believed to be well jus- 
tified for the clusters with closed shells [23|, i.e. for 
TV = 8,18,20,34,40,58,92,..., and in the present work 
we will assume that this condition is satisfied. We un- 
derline, however, that our results are of general char- 



FIG. 1: Metallic cluster formed by alkali atoms Na-8 in mi- 
crocavity. The cavity photons of the frequency ojq undergo 
multiple reemissions and reabsorbtions by a collective exci- 
tation of the cluster forming light- matter hybrid polariton 
eigenmodes of the system. 



acter and should be qualitatively the same for clusters 
with partially filled shells and for clusters consisting of 
the non- alkali atoms. The computation procedure will 
become more tricky in these latter cases [24 1. 

Let us consider an alkali-metal cluster with transition 



frequencies a;, 



— = ~ ^a) where the 

indices a and j correspond to the occupied and unoccu- 
pied single electron states respectively. In the case when 
electron- electron interactions are neglected or treated 
within Hartree- Fock approximation, each single electron 
transition is coupled to the external photons indepen- 
dently. As oscillator strengthes of these individual tran- 
sitions are small, their coupling with cavity mode will 
result in low values of Rabi splitting (of the order of mag- 
nitude of those observed in QD- cavity systems). 

However, the account of the many- body corrections 
can change the picture dramatically. Coulomb interac- 
tion can lead to the formation of the collective excitation 
of the cluster known as giant plazmonic resonance and 
having huge oscillator strength (several orders of magni- 
tude b igger th en oscillator strengthes of individual transi- 
tions) [1 



221 . Consequently, one can expect high values 



of the Rabi splitting in coupled cluster- cavity system. 

Formalism. The interaction of the electromagnetic 
field with material objects is described by dipole matrix 
elements of the transitions daj- The account of many 
electron processes results in renormalization of dipole 
matrix elements which become frequency-dependent and 
are denoted as D^jioj) in our further discussion. 

In diagrammatic representation, the equation for 
Daj{oj) for the system we consider is shown at Figl2l 
and graphical representations of renormalized and bare 
dipole matrix elements of photoabsorbtion D^jioj) and 
daj on Figl3l At the right hand side of Fig 12] the first 
diagram corresponds to the direct photoabsorbtion, the 
second two diagrams accounts for the influence of dy- 
namical polarizability in random phase approximation 
with exchange (RPAE) and give rise to giant plasmon 
resonance, and the last one accounts for the multiple re- 
emissions and re-absorptions of the cavity photon. The 
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FIG. 2: Diagrammatic representation of the equation for the 
renormalized dipole matrix element of the transition. The 
straight lines with arrows correspond to the electrons in the 
cluster at vacant and filled orbitals j, a, the wave lines to the 
Coulomb interaction, the doted lines - to the cavity photons 
of frequency ojo. 
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FIG. 3: Diagrammatic representation of renormalized and 
bare dipole matrix elements of the transition 



straight lines with arrows correspond to the electrons in 
the cluster at vacant and filled orbitals j, a, the wave 
lines to the Coulomb interaction, the dotted lines - to 
the cavity photons of frequency ujq. 

Using the standard rules of the evaluation of the Feyn- 
mann diagrams, one can easily find that renormalized 
dipole matrix elements can be found from the following 
set of non-linear algebraic equations: 



a' J' 

where 

M.,„,,,,(DH) = (-''^'\^o\a,j)-{af\Vc\a',j) ^^^^ 

2a;o^*/^v {Da'j' {^))9aj {Daj {oj)) 

H . 

(cj^ — (jJq -\- 2iujoj) {uj — uJa'j' + ^r) 

In the above formula {a' j'\Vc\oL-, j) denotes the matrix 
element of the Coulomb interaction between the elec- 
trons, the linewidths F and 7 correspond to the finite 
lifetime of the one-electron excitations and cavity pho- 
tons, Qaj are the matrix elements of the interaction be- 
tween the transition a ^ j and cavity mode, which are 
proportional to the renormalized dipole matrix elements 
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of the transitions and can be estimated as [2. 



where V ^ (Ao/2)^ denotes the cavity volume. 

Note, that solving the system of equations (j2j) is equiv- 
alent to accounting the Coulomb RPAE diagrams and 
coupling between transitions in cluster and cavity mode 
up to the infinite order. Together, they describe the 
formation of the hybrid plasmon- photon excitations in 
cluster- cavity system. 

The response of the system (or photoabsorption coeffi- 
cient) is proportional to the imaginary part of the dipole 
dynamical polarizability 
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FIG. 4: Diagrammatic representation of Dyson equation for 
the Green's function of the cavity photon. The double dashed 
line corresponds to the renormalized Green function of the 
photon, the single dot line- to the Green photon of the bare 
photon. The dot on the diagram corresponds to the renormal- 
ized matrix element of the cluster- cavity coupling accounting 
for many- body interactions in the cluster but neglecting mul- 
tiple re-emissions and re-absorbtions. 



which can be found using standart formula 



(5) 
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The energies of the eigenstates of the system can be 
found from the poles of the dressed Green function for 
the cavity photon, represented diagramatically on Fig. HI 
The double dashed line corresponds to the renormalized 
Green function of the photon, the single dot line- to the 
Green photon of the bare photon = 2ujo / (cj^ — cJq + 
2iuJoj). The dot on the diagram corresponds to the renor- 
malized matrix element of the cluster- cavity coupling 
gaj{D^j {(j))^ accounting for many- body interactions in 
the cluster but neglecting multiple re-emissions and re- 
absorbtions of the cavity photon, as latter processes are 
already accounted for in the diagram for in the Dyson 
equation for the dressed photon (FigH]) and should not 
be counted twice. Mathematically, the values of D^^^{uj) 
can be found from the system of the equations analogical 
to those represented in diagrammatic form at FiglSl but 
with last diagram in the right hand side being absent. 
The Green function of the dressed photon reads 



G. 
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where 



^ |ff.,(DgH)|^ 
E(c^) = ^ ^ 



iV 



(8) 



The eigenfrequancies of the hybrid eigenmodes can be 
found from the poles of the expression [7] giving the fol- 
lowing transcendent equation: 



2zc(Jo7 ~ 2uoqT^ {uo) = 



(9) 



Before we proceed with presenting the results of the 
numerical modelling of our system, let us analyze some 
simple limiting cases. 

First, if the cavity is absent, electron- electron inter- 
actions renormalize dipole matrix elements of the tran- 
sitions D^^^iuj) which become frequency dependent and 



form sharp peaks at the frequency uo = uOpi correspond- 
ing to the characteristic frequency of giant plasmon res- 
onance isl-liil (see inset on Fig. 5). 

On the other hand, if one neglects Coulomb interac- 
tions (which corresponds to the retaining of the first and 
last diagrams only in the diagrammatic equation pre- 
sented on Fig|2j) and considers the coupling of the in- 
dividual single transition to the cavity mode, the equa- 
tions for renormalized dipole matrix element and dressed 
photonic Green function can be easily solved and read: 



Do,j{uj) 
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The poles of these expressions determine the energies of 
the new hybrid eigenstates of the system. The condition 
of achievement of the strong coupling regime character- 
ized by the mode and anticrossing at the resonance (cjq = 
ujocj) is given by a standard expression 4^^ > (F— 7)^ |26|. 

Results. Now let us present the results of the numeri- 
cal modelling of the realistic cluster- cavity system. We 
consider a Na-8 cluster embedded in a photonic cavity in 
the position where the electric field of the cavity mode 
reaches its maximum. The energy of the cavity mode ujq 
is supposed to be close to the energy of the giant plas- 
mon resonance uOpi. The energies of the single electron 
states in the cluster were calculated within jellium model 
using Hartree- Fock approximation. The non-radiative 
linewidths of all individual transitions in the cluster were 
taken to be the same and correspond to F = 0.8 meV. 
The nonlinear system of equations (j2j) and[9]for the dipole 
matrix elements and eigenfrequencies of the system were 
solved using iterative procedure of succesive approxima- 
tions. 

Fig El shows the absorbtion of the cavity cluster system 
as a function of frequency of the external laser excita- 
tion for the case when the energy of the giant plasmon is 
tuned in resonance with the energy of the photonic mode, 
uoq = ujpi. If the lifetime of the photons in microcavity 
r is low (7 ^ h/r is high) we can observe one peak on 
the photoabsorption spectrum (see Fig. [5l note, all the 
plots are normalized to unity). However, the situation 
qualitatively changes if lifetime increases. In this case 
the strong coupling regime is established and two peaks 
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FIG. 5: Photoabsorption crossection spectrum - system re- 
sponse on the exciting emittion of light for different hfetimes 
of the photons (in meV): 0.1 (red/sohd), 2 (green/dashed), 4 
(blue/dotted) and 8 (black/dash-dotted). Rabi spUtting Q.r 
decreases at the photons hfetime being decreased and finaUy 
desappears at 7 8 meV. Inset shows spectral characteristic 
of the renormalized dipole matrix element of the cluster Na-8 
(without a cavity). See text for the details. 



FIG. 6: Eigen frequency real part dependence on the mi- 
crocavity mode wavelength for different lifetimes of the pho- 
tons (in meV): 0.1 (red/solid), 2 (green/dashed) and 4 
(blue/dotted). Illustration of the crossover from strong- to 
weak-coupling regime. When 7 8 meV the splitting Q.r = 
and one can see two solutions corresponding to the bare pho- 
ton and plasmon uji = cjq and 0J2 = 00 d ■ 



are observed in photoabsorbtion. The distance between 
the peaks corresponding to the Rabi splitting increases 
with increase of the lifetime of the cavity mode and can 
reach the values of about 20 meV. This is several orders 
of magnitude bigger then the numbers observed in single 
QD- cavity systems 

Fig. [6] shows the dependence of the real parts of the 
eigen frequencies on the cavity energy ujq for different 
linewidths of the cavity (in meV). The Rabi splitting can 
be found as distance between the two branches at the 
anticrossing point corresponding to u = uq and uj = 
ujpi and it varies from meV (7 = 8 meV) to 17 
meV (for 7 = 0.1 meV). The inset illustrates dependence 
of splitting on 7. Note that the results presented on 
Fig. [5] and [6] are in good agreement with each other. It 
should be noted that the increase of the number of the 
atoms in the cluster N leads to the increasing of the Rabi 
splitting. This increase, however is non proportional to 
N- our estimations give maximal values of about 25 meV 
for Na2o cluster (as compare to about 20 meV for Na^ 
cluster). 

Conclusions. In conclusion, we have analyzed the spec- 
trum of photoabsorbtion of an individual Na-8 cluster 
embedded in single mode photonic microcavity. We have 
shown that in the region of the giant plasmonic resonance 
the regime of the strong coupling between plasmon and 
cavity photon can be achieved, which manifests itself by 
formation of the Rabi doublet and mode anticrossing. 
A cluster can be considered as a {)D object with size 
comparable to those of the QDs. However, due to the 
many-body effects it demonstrates Rabi splittings sev- 
eral orders of magnitude bigger then other {)D quantum 
objects. 
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